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Abstract 

We present a calculation of the hadronic vacuum polarization (HVP) tensor within the framework of Dyson- 
Schwinger equations. To this end we use a well-established phenomenological model for the quark-gluon interaction 
with parameters fixed to reproduce hadronic observables. From the HVP tensor we compute both the Adler function 
and the HVP contribution to the anomalous magnetic moment of the muon, a^. We find a^^^ - 6760 x 10"" which 
deviates about two percent from the value extracted from experiment. Additionally, we make comparison with a re- 
cent lattice determination of a^^'' and find good agreement within our approach. We also discuss the implications of 
our result for a corresponding calculation of the hadronic light-by-light scattering contribution to a^. 

Keywords: 



1. Introduction 

One of the most interesting places to search for new 
physics beyond the Standard Model (SM) is the anoma- 
lous magnetic moment of the muon, a^. It is dominated 
by QED eff'ects, however due to the heavy mass of the 
muon it is also sensitive to other corrections. Aside 
from weak interaction contributions which can be eval- 
uated in perturbation theory, one also has to include ef- 
fects from QCD. Since the latter are intrinsically non- 
perturbative at the scales relevant to the calculation, they 
are much harder to include systematically. 

Experimental efforts at Brookhaven National Lab and 
elaborated theoretical eff'orts of the past ten years have 
pinned down a^ to the 10 " level, leading to significant 
deviations between theory [1] and experiment [2, 3]; 
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Figure 1 : The two classifications of corrections to ttie plioton-muon 
vertex function: (a) liadronic vacuum polaiization contribution to a^. 
Tlie vertex is dressed by the vacuum polarization tensor 11^^; (b) the 
hadronic light-by-light scattering contribution to «^ . 



diagram also dominates the theoretical error of a^ from 
the SM. One obtains [4] 
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Experiment: 116 592 089.0(63.0) x 10"" , (1) 
Theory: 116 591 790.0(64.6) x 10"" . (2) 

This 3.3 cr deviation might be seen as a sign for new 
physics, however confirmation requires that the un- 
certainties of both theory and experiment must be 
reduced yet further. The error on the theoretical 
side is dominated by hadronic contributions involving 
non-perturbative QCD dynamics. The leading order 
hadronic contribution is the hadronic vacuum polariza- 
tion insertion (HVP) shown in Fig. 1(a). At present, this 



[6903.0(52.6)- 100.3(1.1)] X 10" 



(3) 



for the leading and subleading contributions, see also 
Refs. [5, 6] for recent updates. The HVP-tensor (H^v) 
involved in the calculation of the leading order result 
can be obtained from experimental input by recourse to 
the optical theorem; such results can then be regarded 
as being model independent. However, note that mod- 
els may be involved in the analysis or extraction of this 
experimental data, especially in the (dominant) low Q^ 
region [7]. 

The diagram that, in the literature, yields the sec- 
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ond largest theoretical error is that of hadronic light- 
by-light scattering (LBL), Fig. 1(b). It is extremely 
difficult to measure and therefore needs to be deter- 
mined from theory alone. There is a long history of 
different approaches to this problem, see Ref. [1] for 
an overview. Recently, we provided a re-evaluation of 
a'lf^ in the framework of Dyson-Schwinger and Bethe- 
Salpeter equations of QCD [8, 9]. Starting with a phe- 
nomenologically successful model for the quark-gluon 
interaction, we determined dynamically the momentum 
dependent quark propagator, the corresponding meson 
Bethe-Salpeter amplitudes and the quark-photon vertex 
and used these as building blocks for our calculation of 
a'lf^. In contrast to previous approaches, we automat- 
ically included effects in the quark-photon interaction 
that are induced by gauge invariance. This can be seen 
as one of the improvements that DSEs have over typ- 
ical effective approaches to QCD. Our results indicate 
that the theoretical value of Eq. (2) may indeed receive 
additional positive contributions that reduce the discrep- 
ancy with experiment. The precise size of these contri- 
butions, however, will only become clear once we re- 
duce the approximations made in [8, 9]. 

While work in this direction is in progress, we find it 
prudent to elucidate upon and justify our approach via 
a calculation of the hadronic vacuum polarization Il^y. 
Although this quantity in principle need not be deter- 
mined from theory for the purposes of a^, it serves as 
an important testing ground for any approach used for 
calculating hadronic contributions to a^ [10, 11, 12, 13, 
14]. In this letter we provide results for the HVP con- 
tribution to the muon anomaly together with the Adler 
function. We employ the same model and philosophy 
as used recently in our calculation of hadronic light- 
by-light scattering [8, 9]. By comparing to the results 
extracted from experiment and to recent lattice calcula- 
tions [15] we will demonstrate that our approach pro- 
vides meaningful and quantitatively reliable results. We 
also believe that our results serve to address and inval- 
idate an argument made by the authors of Ref. [16]. 
There, one-loop radiative corrections to a^^^ and a^^^ 
in a constituent quark model have been invoked to ar- 
gue against large effects from vertex corrections. While 
their calculation is no doubt correct - within the limi- 
tations of using perturbation theory at strong coupling 
scales - the relevance of their results to the case of g- 
2 seems rather limited. This will be discussed in more 
detail below. 

The outline of the letter is as follows. In section 2 we 
will introduce the hadronic vacuum polarization, start- 
ing with its basic definition and its calculation within 
the functional approach. In section 3 we present the 



framework that we employ in this paper, the Dyson- 
Schwinger (DSE) and Bethe-Salpeter (BSE) equations. 
This is followed by our results and a discussion pertain- 
ing fl^^^ and the Adler function in section 4. Finally we 
summarize and discuss the relevance of our results for 
fl^^^ in the concludina sections. 



2. The Hadronic Vacuum Polarization Contribution 

In the following we give the basic definitions con- 
cerning the HVP tensor, the muon anomaly and the 
Adler function. Throughout this work we will employ 
Euclidean space conventions. 

2.1. Basic Definitions 

The hadronic vacuum polarisation tensor Il^y is de- 
fined as that part of the one particle irreducible (IPI) 
photon self energy that is generated by QCD dynamics. 
It can be obtained from the photon Dyson-Schwinger 
equation 
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where D^y is the full photon propagator, Dj,y the bare 
propagator and Z3 is the photon renormalisation con- 
stant. The hadronic tensor n,,v, specified explicitly be- 
low, can also be seen as the IPI-part of the current cor- 
relator 



n^v(?) = J e'^-^j. 
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with /. = / d'''x and the electromagnetic quark current 
;^ given by 



• 2_ 1- ^ 1_ 2_ 1- ^ 

y/i = 2 MT/i" - :^dyf,d - -sy^^s + -cy^,c - -by^,b 
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Here m, d, s, c and b are the respective quark spinors. 
It follows from the Ward Takahashi identity (WTI) 
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that the HVP tensor is transverse: 
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which serves as a definition of the scalar vacuum polar- 
ization Tl{q^). The quantity Tl{q^) is logarithmically di- 
vergent and has to be renormalized. We choose the con- 
dition n(0) = which leaves the definition of the elec- 
tric charge intact. More details concerning our renor- 
maUzation prescription can be found below. 
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Figure 2: The photon vacuum polarization and its expansion in planar 
diagrams. 
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Figure 3: Dyson-Schwinger equation for the quark propagator. Spec- 
ification of the fully-dressed gluon propagator (wiggley line) and 
quark-gluon vertex (grey blob) defines the truncation scheme. 



These are calculated self-consistently within a rainbow- 
ladder approximation to their DSEs, detailed in the next 
section. Note that the diagram in the last line of Fig. 2 
is an exact representation of the hadronic tensor. The 
truncation takes place on the level of the propagator and 
the vertex. 



Once we have obtained the renormalized HVP scalar, 
HrCq'^), the leading hadronic contribution to a^^'' can 
be calculated [10] 
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where m^ is the muon mass and a - e /An is the fine 
structure constant. 

The Adler function D{q) is defined as the logarithmic 
derivative of the polarization scalar 
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The HVP tensor and the Adler function can be obtained 
independently of the model from experiment, exploiting 
dispersion relations (see e.g. [1,4] for details). 

2.2. Expansion in a functional approach 

In a functional approach the vacuum polarisation ten- 
sor is essentially the photon self-energy. For hadronic 
contributions these photons couple to quarks, which in 
turn couple to gluons. It thus contains a resummation of 
an infinity of diagrams. In the spirit of the IjNc count- 
ing we consider only those diagrams which are planar. 
This infinite subset of diagrams is the same as those con- 
sidered in [8, 9]. The resulting expansion is depicted 
graphically in Fig. 2. The first diagram on the right hand 
side shows gluonic corrections that non-perturbatively 
dress the current quark. The second diagram shows 
gluonic corrections to the quark-photon vertex. Both 
classes of diagrams are indicated in the third diagram, 
showing the complexity of the resummation. These are 
finally written in terms of fully-dressed one-particle ir- 
reducible Green's functions (propagators and vertices 
marked by circles) in the second line of the equation. 



3. Framework 

In the following we summarize the calculation 
scheme employed in this paper; more explicit details 
can be found in Ref. [9]. The Dyson-Schwinger equa- 
tions (DSEs) are exact relations amongst the Green's 
function of a given theory. Since they constitute an in- 
finite tower of coupled integral equations a truncation 
has to be employed to provide tractability. For the cal- 
culation of the Adler function and the muon anomaly, 
we need the quark propagator, quark-photon vertex and 
hadronic vacuum polarisation tensor These are ob- 
tained from their respective DSEs, which we detail be- 
low. 

3.1. The Quark DSE 

We begin with the dressed quark propagator S (p), 



s(p)^Zf{p^){i/+M(p^)y 
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which is characterized by the momentum dependent 
quark mass function M{p^) and the wave function 
Zf{p^). These are obtained as a solution of the quark 
DSE given diagrammatically in Fig. 3. On the right 
hand side the inverse bare quark propagator is given by 
S '^(p) - Znii^ + m) with quark renormalization factor 
Z2 and the bare mass m. The quark self-energy contains 
the gluon propagator, given in Landau gauge as 
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with dressing function Z(k~). In addition the dressed 
quark-gluon vertex Tjjip, q) is required. A simple, yet 
phenomenologically successful approximation of the 
quark-gluon interaction has been suggested by Maris 
and Tandy [17]. Here only the leading Dirac structure of 



the vertex is retained r^(A:^) = J^^^'^ik^) and the dress- 
ing of the Yang-Mills (YM) part of the vertex is chosen 
to depend on the gluon momentum k only. The com- 
bination of the gluon- and vertex-dressing functions is 
then modeled as 

Z(;t2)r™(/t2) = ^(JLok^expi-k'/c^) (12) 
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with m, = 0.5 GeV, t = e^ - 1, y„, = 12/(33 - 2Nf), 
Aqcd = 0.234 GeV, w = 0.4 GeV and D = 0.93 GeV^. 
This model interaction assumes the form of the one loop 
running coupling of QCD at momenta k^ » Aq^^ and 
provides enough interaction strength in the infrared for 
dynamical chiral symmetry breaking to occur 

Combining the DSE with the corresponding Bethe- 
Salpeter equation (BSE) one can determine mesonic 
bound state masses and their decay constants. The 
model parameters co and D are then chosen such that the 
physical value of the pion decay constant is reproduced. 
The quark masses have then been fixed by comparison 
with experimental meson masses in the pseudoscalar 
meson sector, cf. the first set in Table 1. Together 
with the self-consistently calculated quark-photon ver- 
tex (see below) electromagnetic properties such as elec- 
tromagnetic form factors and charge radii can be ob- 
tained [18, 19] that are in good agreement with exper- 
iment. This is also true for heavy flavors as discussed 
in Ref. [20]. Especially important for the calculation of 
the HVP tensor is, however, the behavior of the model 
in the vector meson channel. Here, the deviation to ex- 
periment is on the five percent level, as can be seen from 
the first Une of Table 1 . It is therefore not unreasonable 
to expect that the model provides a good description of 
hadronic contributions to II^v up to potential deviations 
of the order of five to ten percent to the experimental 
value. One possibility to investigate the systematic er- 
ror of the model further, is to fix the bare quark masses 
not with pseudoscalar meson masses, but with the vector 
meson sector. The corresponding values are given in the 
second line of Table 1 . Naturally, this is at the expense 
of the pseudoscalar sector, which reacts quadratically to 
a change in the mass parameters, as opposed to the lin- 
ear change of the vector meson sector. Below, we will 
employ both mass parameter sets in our calculation of 
flj^^^ and estimate the model inherent systematic error 
by a comparison of the results. 

3.2. The quark-photon vertex 

The second ingredient necessary for the determi- 
nation of the hadronic tensor through Eq. (2) is the 
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set I 


3.7 


85 


138 


495 


740 


1080 


set II 


11 


72 


240 


477 


770 


1020 



Table 1: Two choices for the hght bai'e quark masses at ^ = 
(19GeV)^ and the resulting meson masses (in MeV) in the pseu- 
doscalar and vector meson sector. For the heavy quarks we always 
take nic = 827 MeV and mt = 3680 MeV which lead to good results 
for charmonia and bottomonia in the pseudoscalar and vector channel. 



fully dressed quark-photon vertex. This quantity is ob- 
tained self-consistently from its inhomogeneous Bethe- 
Salpeter equation 
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where r ^ q - k and ^ ^ J ^ 



We show the BSE 
pictorially in Fig. 4. Again we use the ladder trunca- 
tion ensuring that both the axial-vector and vector Ward 
Takahashi identities (WTIs) are satisfied. To this end 
the quark-gluon interaction in (13) needs to be the same 
as the one in the quark DSE. 

In general, the quark-photon vertex can be decom- 
posed into twelve covariants 
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where V^ are the covariant tensor structures and 
A^'\P,k) are non-trivial scalar dressing functions that 
contain the non-perturbative dynamics. The photon mo- 
mentum is P, with k the relative quark momentum such 
that the incoming and outgoing quark momenta are k-t - 
k + Pjl. Ball and Chiu suggested to separate the vertex 
into the transverse parts V^ ' ' ' with P^V^, '' ' - and 
four non-transverse components V^ ' " . The latter ones 
are completely fixed in terms of the quark dressing func- 
tions M{p-) and Zf{p-) by the WTIs and the demand of 
regularity [21]. This part of the quark-photon vertex is 
also called the Ball-Chiu- or BC-vertex. 
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Figure 4: Inhomogeneous BSE for the quark-photon vertex. 



The additional eight components of the trans- 
verse part are determined numerically through a self- 
consistent solution of the BSE, Fig. 4. It is an important 
property of this equation that it generates dynamically 
vector meson bound-state poles. The idea of vector me- 
son dominance then corresponds to the suggestion that 
the transverse part of the quark-photon vertex provides 
the leading contributions in a calculation at hand. While 
this is indeed correct for some observables, as for ex- 
ample for a^^'' below, other examples are known where 
this idea is not correct and sizable contributions from 
the BC-part of the vertex occur, see e.g. [22]. 

3.3. The hadronic vacuum polarization 

Finally we give some details regarding the calcula- 
tion of the hadronic vacuum polarization. Within the 
truncation scheme proposed above, the hadronic tensor 
is given by 
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where q± - q ± P/2 and Z2 is the quark wave function 
renormalisation. The scalar function n(P^) is obtained 
via Eq. (7). This quantity is logarithmically divergent 
and so requires renormalisation. We apply the condition 
nR(0) - through the subtraction 



nR(p2) := n(p2) - n(0) 



(16) 



which effectively amounts to adjusting the constant Z3 
in Eq. (4) appropriately. In addition we need to take 
care of quadratic divergences that appear through our 
use of a hard numerical cutoff. These can be subtracted 
at /?^ = or projected out using the method of Brown 
and Pennington [23]. Both procedures agree very well. 

To check our numerics, we first evaluated the per- 
turbative QED one-loop result (see e.g [24]) and found 
excellent agreement. In particular we checked that the 
calculation of IIr was independent of the cut-off. As a 
further check, we evaluated the electron loop contribu- 
tion to fl^ via Eq. (8) by replacing the propagator and 
vertices with their tree-level values. We reproduced the 
well known result 0™'=?°'-'=-!°°? ~ 5.994 x 10"'' [1] on 
the sub per mille level. For our general calculations 
with dressed momentum dependent quark propagator 
and quark-photon vertex we estimate a numerical er- 
ror of roughly two to three percent due to the uncer- 
tainties related with the renormalisation procedure dis- 
cussed above. 

Below we present the results of our calculation for 
the Adler function as well as the anomalous magnetic 
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Figure 5: The quark mass functions of the 11/ d, s, c and b quarks 
obtained from the quark DSE. The dashed vertical line represents the 
renormalization point//- = (l9GeV)^. 



moment of the muon a^. We use the Maris-Tandy inter- 
action with the two different parameter sets discussed 
above. We solve the quark DSE, Fig. 3, for five quark 
flavors M, d, s, c and h, and work in the isospin symmet- 
ric limit m„ = mj. The resulting quark mass functions 
are shown in Fig. 5. They dynamically connect the in- 
frared constituent quark mass region with the ultravio- 
let current quark mass region and thus provide a unified 
approach to both pictures. Note that our quark agrees 
qualitatively with lattice calculations [25]. 

Once the quarks are obtained we solve for the quark- 
photon vertex, Eq. (13). Here, no additional approxima- 
tions are made, i.e. we take into account all twelve ten- 
sor structures and the full momentum dependence of the 
vertex. This is done for each flavor separately and hence 
we can calculate H^v using Eq. (15) which sums over 
all quark flavors. With n^,v at hand we can obtain the 
hadronic contribution to the anomalous magnetic mo- 
ment of the muon via Eq. (8), and the Adler function 
fromEq. (9). 



4. Results 

In Fig. 6 we show our result for the Adler function as 
calculated using parameter set II of Table 1, together 
with the result from dispersion relations [1, 4]. The 
Dyson-Schwinger solution describes the data very well 
in the non-perturbative region Q < \ GeV. We also see 
that in the asymptotic ultraviolet limit the solution fol- 
lows the result from the dispersion relations. The differ- 
ences between set I (not shown) are limited to the slope 
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Figure 6: The Adler function obtained from DSE's for the Maris- 
Tandy model together with the dispersion relation results from [26, 4]. 
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Figure 7: The Adler function obtained in the MT model defined in Eq. 
(12) via Eq. (15) with diiferent vertex dressings. 



of the function in the low momentum region, which is 
most sensitive to the mass of the vector meson (see fig. 
7). Note in addition that most of the contributions to 
a^^^ come from the region around the muon mass and 
that the integration of Eq. (8) saturates between 0.5 and 
1 GeV. From the Adler function we therefore expect 
similar results for a^^'' for both parameter sets with 
small deviations on the level of ten percent. 

Before we discuss our results for a^^'' we take a 
closer look at the impact of the transverse parts of the 
quark-photon vertex as compared to its non-transverse 
Ball-Chiu (BC) structure. In Fig. 7 we compare the full 
results with the one using the BC-part alone or even ne- 
glecting all vertex dressing altogether Clearly, the bare 
and Ball-Chiu vertices do not provide sufficient contri- 
butions to the Adler function, yielding functions that are 
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bare 


BC 


transverse 


full 


set: 


760 


1280 


6160 


7440 


set II 


720 


1120 


5640 


6760 



Table 2: The leading order HVP contribution to a^ as obtained by our 
two sets of bare quark masses for different truncations of the quark- 
photon vertex. 



only half the height of the full vertex result in the in- 
frared. Only the full vertex calculated from its inhomo- 
geneous Bethe-Salpeter equation contains vector meson 
poles dynamically in its transverse structure. Obviously 
these are essential to describe the data correctly. 

This sensitivity to the vector meson sector is espe- 
cially seen in a^^''. For the two mass parameter sets I, 
II of our model and the full quark-photon vertex we find 
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As expected, our first mass parameter set yields a value 
for a^^^ which is too large by about eight percent, due 
to the fact that our vector meson for this parameter set 
is slightly too light and can thus be excited from the 
vaccum too easily. This, however, is already a reas- 
suringly good result for a calculation performed with 
standard parameters without adjustment. Changing our 
input mass parameters to values that are matched to the 
vector meson sector improves our value for a^^^ such 
that deviations with experiment fall below three percent. 
We regard this agreement as a clear signal that our ap- 
proach accurately contains the physics relevant for the 
hadronic contributions to a^, which entails that indeed 
the dynamics associated with the vector meson pole, to- 
gether with gauge invariance, are the two essential in- 
gredients. 

Next we examine the dependence of a^^'' on the 
quark-photon vertex used in Eq. (15). The results can be 
found in Table 2. As expected from our results for the 
Adler function, most of the contribution to a^^'' comes 
from the transverse parts of the vertex containing the 
vector meson poles. Here also most of the differences 
between our parameter sets I and II occur. However, 
there are also sizable contributions from the gauge or 
Ball-Chiu part of the vertex and only the use of the full 
vertex gives satisfying results for a^^''. Once more, this 
emphasizes the interplay of contributions related to res- 
onances and those demanded by gauge symmetry. 

Finally we look closer at the dependence of a^^'' 
on the quark mass. This behavior is conveniently 
parametrized by plotting against a scheme independent, 
physical mass such as for example the pseudoscalar or 
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Figure 8: The mass dependence of a" X 10 , for two flavours, 
plotted wrt the mass of the light vector meson. Shown is the data 
from this work (DSE) together with a fit (a^ oc My " ). In addition 
we show recent data from the ETMC lattice collaboration [15]. 



vector meson mass my- Both of these can be determined 
in our approach via a solution of their corresponding 
Bethe-Salpeter equations. In Fig. (8) we show our re- 
sults for a^^''{mv) compared to a recent lattice study 
of the ETMC-collaboration [15]. Overall we find very 
good agreement between the two approaches, with our 
values inside their error bars. The same level of agree- 
ment is seen between our calculation and the new lattice 
determination for A^y = 2 H- 1 flavour QCD presented in 
Ref. [27]. 

5. Discussion 

Our results for a^''^ clearly show the importance of 
dressing effects in the quark-photon vertex. Here, par- 
ticularly relevant are its dynamically generated vector 
meson poles in the transverse part of the vertex. How- 
ever, we wish to emphasize again that this importance 
crucially depends on the kinematic and dynamical de- 
tails of the problem at hand. For example, the trans- 
verse parts of the vertex contribute towards only half 
of the pion charge radius [18], whilst in the pion pole 
approximation of the light by light contributions to g-2 
they constitute only a thirty percent effect as compared 
totheBCpart [8, 9]. 

It is thus very dangerous to transport expectations 
based on one process blindly to another; explicit cal- 
culations should always be preferred. In this work, we 
have performed such a calculation for a^^'' by including 
both the BC- and transverse parts of the vertex explic- 
itly. For a^^^ in Refs. [8, 9], the algebraic complexity 



forced us to consider initially only the BC part of the 
vertex, with transverse parts estimated from other model 
calculations. Preliminary results for aff^ with the full 
vertex have been presented at [28], and show that gauge 
effects still dominate. The details of this will be pre- 
sented in a future work. 

We also wish to discuss the arguments made in 
Ref. [16]. There, a constituent quark model with mo- 
mentum independent masses has been combined with 
a perturbative evaluation of gluonic corrections. Cor- 
responding results for a^^'' and a^*^ have been com- 
pared. The authors point out that neglecting radiative 
corrections, they need unphysically small constituent 
quarks masses to reproduce the experimental value for 
a^^''. From our results we can see clearly that this is 
merely the result of compensating for dynamics that are 
absent in the quark-photon interaction of that model. 
The authors take note of that fact and argue that this 
very light constituent quark mass effectively includes 
the -y - p-coupling gp via Mg oc Mp/gp. This simple 
relation might however be inappropriate for very dis- 
similar kinematics. In addition the authors of Ref. [16] 
find very large corrections when they include radiative 
corrections on the one-loop level. They observed that 
these corrections could be absorbed into a change of 
the constituent quark mass with stable results for a^^'' 
and Op^^. Based on this result the authors suggest that 
dressing effects in the quark-photon vertex of the full 
theory should be small. We disagree with this conclu- 
sion. First of all, it is dangerous to interpret a truncated 
perturbative expansion that features both a large expan- 
sion parameter and large expansion coeflicients. Sec- 
ond, non-perturbative features such as the formation of 
bound-states (as generated dynamically by the vertex) 
are absent in their calculation. Thus we take their re- 
sults as a hint that (infinitely many) higher contributions 
are important and should be included as consistently as 
possible, since any finite order pQCD cannot give satis- 
factory answers. 

We have done exactly this in our calculation. As a re- 
sult we found that the leading order contribution comes 
from vector meson (VM) poles accounting for roughly 
80 % of Op^'', with the remainder coming from correc- 
tions induced by gauge invariance. The dynamics of the 
VM poles are thus important but not the whole story. 
This tells us that an effective model that features only 
VM exchange should be a good approximation, but will 
miss out on other important contributions that cannot be 
integrated by reshuffling of contributions. Similarly, a 
constituent quark loop approach would not contain any 
dynamical degrees of freedom relating to vector meson 



exchange. This is in contradiction to what is observed 
both on the lattice and in our Dyson-Schwinger calcu- 
lation, and thus the constituent quark model cannot be a 
satisfactory description of the process at hand. 

Finally, we believe that the good agreement of our 
results for a^^'' with experiment and with lattice cal- 
culations adds credit to our corresponding approach to 

^,LBL 



6. Summary 

We calculated the hadronic vacuum polarization us- 
ing the method of Dyson-Schwinger equations, taking 
into account the five lightest quark flavors. As input 
we used a phenomenologically successful model for 
the quark-gluon interaction together with the rainbow- 
ladder truncation. The parameters of these interactions 
as well as the quark masses were fixed by meson ob- 
servables such as masses and decay constants, without 
additional fine-tuning. We determined the quark-photon 
vertex from its inhomogeneous Bethe-Salpeter equa- 
tion in the same approximation and subsequently cal- 
culated the hadronic vacuum polarization tensor From 
these we obtained results for the anomalous magnetic 
moment of the muon a^ ' as well as for the Adler 
function. Both quantities agree well with model inde- 
pendent results extracted from experiment. In particu- 
lar, the Adler function is reproduced very well in the 
strictly non-perturbative region at small momenta. We 
have shown that one requires a description in terms of 
dynamical quarks interacting through non-perturbative 
gluons in order to achieve this level of accuracy. 

Consequently we find results for the muon anomaly 
in good agreement with other determinations. Our best 
result using the quark mass parameter set II is 



flHVP,LO ^6 760x10- 



(19) 



This can be compared to the leading order result quoted 
in Eq. (3), 6 903.0(52.6) x lO"". The difi^erence is at 
the level of two percent. A comparison with the result 
^HVP,LO _ ^44Q ^ ^Q-ii Q|3t^jjjg(j ^jfjj Qjjj parameter 

set I may serve as an estimate for the systematic uncer- 
tainty of our model of roughly ten percent. We believe 
our approach to the hadronic light-by-light scattering 
contribution [8, 9], which employs the same truncation 
scheme, will ultimately lead to results with similar pre- 
cision. However, note that in Ref. [8, 9] the full quark- 
photon vertex was not yet included in the quark-loop 
due to its algebraic complexity. Improvements along 
this direction are underway. 



7. Acknowledgments 

We thank A. E. Dorokhov and A. E. Radzhabov for 
helpful discussions. This work was supported by the 
DFG under grant No. Fi 970/8-1, by the Helmholtz- 
University Young Investigator Grant No. VH-NG-332 
and by the Helmholtz International Center for FAIR 
within the LOEWE program of the State of Hesse. RW 
would also like to acknowledge support by the Austrian 
Science Fund FWF under Project No. P20592-N16, 
and by Ministerio de Educacion (Spain): Programa Na- 
cional de Movilidad de Recursos Humanos del. Plan 
Nacional de I-D-Hi 2008-201 1. 



References 



[1 

[2 

[3 

[4: 
[5 

[6: 

[7 

[8 

[9 

[10 

[11 

[12: 

[13 

[14: 

[15 

[16: 

[17 
[18 

[19 
[20 
[21 
[22 
[23 
[24" 
[25 

[26: 
[27 



F. Jegerlehner and A. Nyffeler, Phys. Rept. 477, 1 (2009) 
[arXiv:0902.3360 [hep-ph]]. 

G. W. Bennett et al. [Muon G-2 Collaboration], Phys. Rev. D 73 
(2006) 072003 [arXiv:hep-ex/0602035]. 

B. L. Roberts, Chin. Phys. C34 (2010) 741-744. 
[arXiv:1001.2898 [hep-ex]]. 

F. Jegerlehner, Berlin, Germany: Springer (2008) 426 p 

M. Davier, A. Hoecker, B. Malaescu et al, Eur. Phys. J. C71 

(2011) 1515. [arXiv:1010.4180 [hep-ph]]. 

K. Hagiwai'a, R. Liao, A. D. Martin, D. Nomura, T. Teubner, 

[arXiv:l 105.3149 [hep-ph]]. 

M. Benayoun, P. David, L. DelBuono, F. Jegerlehner, 

[arXiv:l 106.1315 [hep-ph]]. 

C. S. Fischer, T. Goecke, R. Williams, Eur. Phys. J. A47 (201 1) 
28. [arXiv: 1009.5297 [hep-ph]]. 

T. Goecke, C. S. Fischer, R. Williams, Phys. Rev. D83 (2011) 

094006. [ai-Xiv: 1012.3886 [hep-ph]]. 

E. de Rafael, Phys. Lett. B 322 (1994) 239 [arXiv:hep- 

ph/9311316]. 

E. Pallante, Phys. Lett. B341 (1994) 221-227. [hep- 

ph/9408231]. 

J. S. Bell, E. de Rafael, Nucl. Phys. Bll (1969) 611-620. 

A. E. Dorokhov, Phys. Rev. D70 (2004) 094011. [hep- 

ph/0405153]. 

M. Delia Morte, B. Jager, A. Juttner, H. Wittig, AIP Conf Proc. 

1343 (2011)337-339. [arXiv:1011.5793 [hep-lat]]. 

X. Feng, K. Jansen, M. Petschlies, D. B. Renner, 

[arXiv:1103.4818 [hep-lat]]. 

R. Boughezal, K. Melnikov, [arXiv:1104.4510 [hep-ph]]. 

R Maris and P C. Tandy, Phys. Rev. C 60 (1999) 055214. 

P Maris and P C. Tandy, Phys. Rev. C 61, 045202 (2000) 

[arXiv:nucl-th/9910033]. 

M. S. Bhagwat and P Maris, Phys. Rev. C 77, 025203 (2008). 

P Maris, AIP Conf Proc. 892, 65 (2007). 

J. S. Ball and T W. Chiu, Phys. Rev. D 22, 2542 (1980). 

P Maris and P C. Tandy, Nucl. Phys. A 663, 401 (2000). 

N. Brown, M. R. Pennington, Phys. Rev. D39 (1989) 2723. 

V. P Nair, New York, USA: Springer (2005) 557 p" 

P. O. Bowman, U. M. Heller, D. B. Leinweber, M. B. Pai'appilly, 

A. G. Williams and J. b. Zhang, Phys. Rev. D 71 (2005) 054507 

[arXiv:hep-lal/0501019]. 

S. Eidelman, F. Jegerlehner, A. L. Kataev and O. Veretin, Phys. 

Lett. B 454, 369(1999). 

P. Boyle, L. Del Debbio, E. Kerrane and J. Zanotti, 

arXiv:l 107.1497 [hep-lat]. 



[28] R. Williams, talk at the Workshop on the hadronic light-by-light 
scatteiing contribution to the muon magnetic anomaly, INT, 
Seattle 02/28-03/04 2011. 



